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We present a comparative numerical study of the ordered and the random two-dimensional 
sine-Gordon models on a lattice. We analytically compute the main features of the expected 
high temperature phase of both models, described by the Edwards- Wilkinson equation. We 
then use those results to locate the transition temperatures of both models in our Langevin 
dynamics simulations. We show that our results reconcile previous contradictory numerical 
works concerning the superroughening transition in the random sine-Gordon model. We also 
find evidence supporting the existence of two different low temperature phases for the disordered 
model. We discuss our results in view of the different analytical predictions available and 
comment on the nature of these two putative phases. 
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I. INTRODUCTION 
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The location and characteristics of phase transitions 
in disordered media constitute a long-standing and con- 
troversial question, particularly in more than one spatial 
dimension The question becomes even more diffi- 
cult if the disorder is not very weak; then, new, non- 
trivial behavior is commonly found, involving features 
such as aging, ergodicity breaking, extremely slow dy- 
namics, complicated energy landscapes, etc.; major ex- 
amples of this are spin glasses and structural glasses 
10, ||]. In this context, the properties of crystal surfaces 
growing on disordered substrates, frequently described 
by a two-dimensional random-phase sine-Gordon model 
(RsGM), have attracted a lot of attention in the last 
decade 
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Interestingly, the same model describes 
many other relevant physical problems, such as randomly 
pinned flux lines confined in a plane [|[ |2^, p5| , 
vortex lines in planar Josephson junctions [pq , charge 
density waves [ p7[ , and commensurate-incommensurate 
transitions ||2^ . 

In spite of those efforts, the phase diagram and main 
features of the RsGM are not clear yet. To summarize 
what is known, we refer for comparison to the ordered 
sine-Gordon model (OsGM), which is rather well under- 
stood (see, e.g., ||l|, |2^, |3^, |3l|, |3^, |3^ and references 
therein). The hamiltonian for the OsGM and the RsGM 
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where n.n. stands for the nearest neighbors of site r. 
The OsGM corresponds to /i*^°)(r) = 0, whereas the 
RsGM is defined by choosing the quenched disorder vari- 
ables h^'^\r) randomly from a uniform distribution in 
[0, 27r]. We will discuss our work in terms of growth 
on fiat (OsGM) or disordered (RsGM) substrates (see 
p^ , p5| , [2^ , p7| for other physical interpretations): Ac- 
cordingly, h{r) is a continuous variable representing the 
height of the growing surface at site r of the lattice, and 
the cosine term is a potential energy making integer (i.e., 
multiples of the crystalline lattice constant) values of 
the height energetically favorable. We will consider two- 
dimensional (2D) square lattices, i.e., r = (ri,r2), with 
N — L X L sites. As first shown by Chui and Weeks Q 
(see also paj ) by means of a renormalization group (RG) 
approach |36[, the OsGM posseses a Kosterlitz-Thouless 
P7| type topological transition between a low tempera- 
ture, fiat phase and a high temperature, rough phase, the 
latter being described by the Edwards- Wilkinson (EW) 
equation pS] , i.e., the diffusion equation with additive 
white noise (see below). Above the so called roughening 
temperature (T/{), thermal fiuctuations effectively sup- 
press the effect of the cosine potential, and the surface 
becomes free, described only by the kinetic part of the 
hamiltonian (^. As will be shown below, the most im- 
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portant measurable consequence of this is that the width 
of the interface, 
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scales (in 2D) as ^ In L in the asymptotic regime. 

In contrast with the clear picture for the OsGM, there 
are very few generally accepted results for the RsGM. 
One of them is that there must be a roughening tem- 
perature above which the potential effectively vanishes 
(much as in the case of the OsGM) leading to EW be- 
havior. Apart from this, theoretical predictions about 
the low temperature phase largely disagree (a good sum- 
mary is given in the third paper in Ref. §): While 
RG calculations predict a superrough low temperature 
phase, with up' ~ In^ L, replica-symmetry breaking vari- 
ational approaches lead to '--^ InL independently of 
temperature. Numerical simulations were not very con- 
clusive either: Batrouni and Hwa ||^ did not find ev- 
idence for an equilibrium phase transition in Langevin 
dynamics, although Monte Carlo simulations by Rieger 
[|6[ showed a transition from a superrough phase to a 
EW phase for stronger potentials [larger Vq in Eq. (|^)] 
than those used by Batrouni and Hwa. Subsequent nu- 
merical work n, [l2[ ^ |l^ presented more evidence of 
superrough (In Ljbehavior, albeit with large quantita- 
tive discrepancies with the predictions of RG theories. 
Finally, a number of works using special optimization al- 
gorithms |l^ or direct numerical simulations 
strongly supported superrough behavior at zero temper- 
ature. Very recently ||2^, simulations of a related model 
provided more evidence of In^ behavior at finite temper- 
atures, although this model did not allow study of the 
transition. In summary, most researchers believe that 
there is indeed a superrough low temperature phase in 
the RsGM, but its nature (glassy or not), the transition 
temperature, and its dependence on the model parame- 
ters remain unclear. 

In this paper, we attempt to shed light on these issues 
by simultaneously studying the OsGM and the RsGM 
in different regions of their phase diagram. As we will 
show below, it turns out that the potential strength, Vb, 
crucially determines the model features. In addition, it 
is also natural to ask about the intensity of the disorder. 
How does the model phenomenology change if the disor- 
der takes values in [0, e] with e < 2tt1 The importance of 
these points can be clearly seen in the Langevin equation. 
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where ri{r,t) is a gaussian white noise of zero mean and 
correlations obeying the fluctuation-dissipation theorem, 
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where T stands for the temperature of the system, (...) 
indicates thermal averages (over rf) and Boltzmann's con- 
stant is set to fcs = 1. For our hamiltonian, Eq. dl^), the 



Langevin equation is 
dh{r,t) 



dt 



=V'h{r,t) - 
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(5) 



and changing variables according to h{r,t) — h{r,t) — 
h'^^^(r) (i.e., the height referred to the substrate), we find 
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where F{r) ~ \/^h^-^\r)/e. In this form, the disorder 
is a random (correlated in space) chemical potential act- 
ing on a surface growing on a flat substrate. If we think 
of the roughening transition for the OsGM in terms of 
the interplay between the temperature T and the energy 
scale introduced by the potential, Vq, when e 7^ 0, we 
have another energy scale, e, which can modify the uni- 
versal features of the roughening transition or even give 
rise to novel thermodynamical transitions. 

Having the above issues in mind, we discuss our work 
according to the following scheme: Section II presents 
an analytical study of the EW equation and other sta- 
tistical mechanics results about the energy and rough- 
ness. Section HI deals with the main part of our work, 
namely Langevin dynamics simulations of the OsGM and 
the RsGM, beginning with Vq — I and e — 2tt and sub- 
sequently analyzing the model behavior for different Vq 
and e. Finally, in Sec. IV we present and discuss our 
conclusions and indicate future directions. 



II. ANALYTICAL RESULTS 

A Linear theory: The Edwards- Wilkinson equation 

According to the RG approach the high temper- 
ature phase of the OsGM obeys the EW equation pm , 



dh{r,t) 
dt 



V^h{r,t) +r){r,t). 



(7) 



Equation (Q) can be solved by means of the Fourier de- 
composition (see pa] for details). 



(8) 



where q = ^k, ki ^ 0, . . . ,L — lis the reciprocal vector. 
The structure factor can then be shown to be 

S{q) = {h,h^g)^T^^^^^, (9) 
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ujq being the 2D EW discrete dispersion relation 



ujq = 4sin2 



4 sin' 
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From S{q) we can obtain the relevant magnitudes, such 
as the total roughness 
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the correlation function, 

C{r) = /^J2^h{s + r)-h{s)f 



the total slope, 

and the slope difference correlation function, 

4 ^ 
= 72 E -^(^^ Et-^ ~ cos(g,ri)][l - cos(9. 



(12) 
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From Eq. (|) we can find an estimate for the time 
needed to reach saturation, tx, and the dynamic expo- 
nent, z: We compute the time that the structure factor 
needs to be within a 1% of its saturated form for the slow- 
est Fourier mode, that with the lowest l^j, \q\ — 2n/L: 

tx~3xlO~2L2, (15) 

implying z — 2. For the saturated roughness, we obtain 
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which cannot be computed exactly but, for large L, can 
be approximated changing the sum by an integral and 
the sine functions by their arguments, arriving at 
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yielding a roughness exponent a 0. As for the total 
slope, it tends to a value independent of L, 



T, 



(18) 



whereas for the correlation functions, we find for large r 



T 

C(r) ~ — Inr, 



(19) 



T 



G(r) ~ 2s + In ^^^^ ^ _ 4rl][{r^ + 1)2 - Ar^] 
~ 2T. (20) 

B Other results for the energy and the roughness 

At equilibrium, the partition function for the OsGM is 



Z = 



Y[dh{i 



(21) 



where fi = T ^. Expanding for high temperatures | 
means rewriting Eq. (E^) as a series in powers of /3Vo: 
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where Zq is the partition function of Hq, the free hamilto- 
nian [i.e., Eq. (n|) without the potential term]. By means 
of this expansion we obtain 
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X ((^cos/i(r))"cos/i(r'))ifo> 



(23) 



where subdominant terms such as /3" exp(— A//?) have 
been neglected, and only terms of the order /3" have been 
kept. The expression above can be put in the form 

{cos h{r'))H = 2^ ^-:ttt7T-^ I X 
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(/31^o)'" 



(24) 
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(17) yielding for the approximate energy per site 
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For the roughness, we have to compute 
mr) - hf)H - f X 



(26) 
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assuming that at equihbrium h = Q. Neglecting again 
subdominant terms, we find 



(27) 



and hence 
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Finally, at low temperatures the height exhibits only 
small deviations from h — Q, and therefore we can ap- 
proximate the hamiltonian as 



i 



(29) 



i.e., 2N quadratic terms, each one of which, according to 
the equipartition theorem |^^, contributes with T/2 to 
the energy value. Taking into account the global factor 
1/2, we conclude that at low temperatures the energy of 
the OsGM is approximately E = T/2. 



III. LANGEVIN DYNAMICS RESULTS 

A Numerical simulation details 

We have integrated the Langevin equation (^) corre- 
sponding to the hamiltonian (|l|) , on L x L square lattices 
with periodic boundary conditions, using a stochastic 
second order Runge-Kutta method in some cases, 
we have repeated the simulations with a Heun method 
[Q, with excellent agreement between both procedures. 
We are therefore sure that our results are not an artifact 
of our numerical method, a conclusion further reinforced 
by the agreement with the theoretical expectations of 
the previous section as we discuss below. The simula- 
tions reported in this paper were carried out with a time 
step At — 0.01 on lattices of sizes L = 64, 128 and 256. 
It is important to stress that we did not perform aver- 
ages over the quenched noise in the RsGM; however, we 
checked that the outcome of our simulations did not de- 
pend strongly on the realization of the quenched noise 
or the initial conditions (flat, h{r) = 0; as the substrate. 



h{r) = h^°\r), or random) by repeating several times a 
number of our simulations. In all cases, simulations con- 
sisted of an equilibration time and a measuring period. 
Eq. (|l^) predicts that the time needed for equilibration 
is tx ^ 500 for L = 128 and ~ 2000 for L = 256, 
and hence we used equilibration times of 5000 and 10000 
units, respectively; afterwards, we let the system evolve 
for an equal period, over which we performed thermal av- 
erages. Equilibration was ensured in all cases by verifying 
that the fluctuations of the energy were gaussian and by 
checking the equality of the specific heat computed from 
those fluctuations and from derivatives of the mean en- 
ergy ||l| , as well as by monitoring the evolution in time of 
the quantities of interest toward a stationary state. As 
an additional test, we compared the imposed simulation 
temperature, arising from the noise term, to that mea- 
sured during the evolution according to the equipartition 
theorem ; both quantities were always found to agree 
within a 0.1%. Finally, we did a few very longer runs, 
whose outcome agreed with that of the shorter runs. 



B Standard RsGM: Vo = 1 and e = 2n 



We begin by discussing our results for the "canonical" 
version of both the OsGM and the RsGM, i.e., hamilto- 
nian (n|) with Vb = 1, as studied (for the ordered case) 
in [p9L^Q]. In those works, the roughening temperature 
was determined by a direct comparison to RG predic- 
tions, looking for the temperature at which the height 
difference correlation function reached a universal (in the 
RG framework) value, with the result that for the OsGM 
Tr « 25 in our dimensionless units. Remarkably, this is 
the RG value for Tfj, which makes very tempting to claim 
that this method indeed yields correctly. Monte Carlo 
simulations of the discrete gaussian model (hamiltonian 
(|l|) with Vb — and h{r) restricted to integer values) by 
Shugard et al. with the same criterion for locating 
the transition yielded similar results. However, as RG 
calculations are perturbative in Vq and carried out on 
the continuum Langevin equation ||2^, ^ ^ , 

it is not obvious that they apply to a discrete model with 
Vb = 1, i.e., of the same order as the kinetic term. In view 
of this, we decided to include the OsGM in this study, 
both to analyze in detail whether the comparison to the 
universal RG prediction for the factor is a good tool to 
find Tn and to compare its high temperature phase with 
that of the RsGM, which should also be of EW type. 

The first quantity we discuss, shown in Fig. |l|, is the 
mean energy of both models. As we see, the results are 
largely independent of the system size, and hence it is 
unlikely that they are affected by finite size effects. The 
plot shows that the mean energy of the OsGM reaches 
the high temperature approximation at Tq = 16 ± 1; on 
the other hand, the mean energy for the RsGM is never 
too far from it, although for temperatures lower than 
Ti = 4 ± 1 the numerical values lie slightly below the 
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high temperature result. At temperatures higher than 
Tq the energies of both models coincide within the ac- 
curacy of our simulations. These results suggest that Tq 
could be the roughening temperature, T^, for the OsGM 
and Ti the superroughening temperature, Tsr, for the 
RsGM, because the EW behavior of the mean energy 
of both models indicates the effective suppresion of the 
sine term by temperature. The inset in Fig. |l| presents 
the specific heat, C^, of both models, exhibiting a well 
defined peak in C^, for the OsGM with its maximum at 
temperature T — 9, much lower than Tq. In contrast, 
we do not observe any peak for the RsGM; there might 
be a peak at T = 3, but the evidence is not conclusive. 
Concerning the peak for the OsGM, we stress again the 
absence of any finite size effect, consistent with a KT 
type transition. Furthermore, we believe that this peak 
is a (Schottky) anomaly similar to that observed in 
2D XY and easy-plane Heisenberg spin models |Q above 
the KT transition. Recall that when mapping the OsGM 
to the XY model the temperature of the former maps 
to the inverse temperature of the latter |Q, hence the 
observation of the anomaly below the possible transition 
temperature Tq. This reinforces our interpretation of Tq 
as the roughening temperature of the OsGM. 

The total roughness of both models, shown in Fig. |^, 
behaves similarly for temperatures higher than Tg, de- 
pending linearly on temperature. In both cases, we see 
that the slope of the roughness depends on the system 
size, as predicted by Eq. (p^). Whereas the approxima- 
tion in Eq. (0) expression yields slopes 0.66 and 0.77 for 
L = 64 and 128, respectively, if we numerically compute 
the exact result, Eq. (|ll|), the slopes turn out to be 0.71 
y 0.82, in excellent agreement with the results of our sim- 
ulations, 0.71 and 0.83. Below Tq, the roughness for the 
OsGM is independent of the size and depends nonlinearly 
on temperature, whereas above Tq we find linear depen- 
dence on temperature and clear finite size effects. For the 
RsGM, the linear behavior extends all the way down to 
Ti, and below Ti the behavior becomes nonlinear. The 
slope of the linear region is approximately the same in 
the first part, from Ti to Tq, and in the second, above To, 
i.e., the whole linear region is well described by the EW 
model. This means that above Tq the linear model de- 
scribes accurately the behavior of the OsGM, and hence 
from now on we identify Tq with the roughening temper- 
ature Tfl, whereas for the RsGM, the same is true of Ti 
and Tsr. Figure || presents also results for the roughness 
susceptibility, Xw, defined as xt, = [((^«^)^) - (w^)^]/^^- 
For the OsGM, Xw exhibits a very clear peak at Tr, and 
above Tr it is the same as for the RsGM; however, this 
magnitude is very noisy and these results must be taken 
with caution. In fact, one could identify a peak for the 
RsGM at Tsr, but different realizations lead to different 
results, in contrast with the peak for the OsGM, which 
is the same for all realizations. Below Tr, the values of 
Xw for the OsGM are independent of the system size, 
whereas above Tr they increase with size without any 
definite scaling. 



Figure ^ depicts the height-difference correlation func- 
tion for the two studied models, and shows that above 
Tr = To and Tsr^ = Ti they behave as predicted by 
the linear theory: The slope of the numerical height- 
difference correlation function is 0.32, indistinguishable 
from the predicted 0.318 by Eq. ([l9|). In addition, the 
correlation functions for the OsGM and the RsGM coin- 
cide, as shown in the plot for the RsGM. We see that the 
behavior of the correlation functions is in full agreement 
with our claims regarding Tr and Tsr, and this is fur- 
ther confirmed by the plot of the slope correlation func- 
tion in Fig. ^ It is important to note that, below Tsr, 
the behavior of the height difference correlation function 
for the RsGM is approximately a squared logarithm, as 
predicted by RG calculations. We postpone discussion of 
this point to the next section. 

Finally we studied another magnitude, namely 



mi = (cos \h{r) - h^"\r)\) 



(30) 



For the OsGM this is the average computed in the pre- 
ceding section, whereas for the RsGM it is the average of 
the cosine of the height referred to the substrate. Figure ^ 
shows our results: The high temperature approximation, 
Eq. ( p^ ) , agrees very accurately with the simulations for 
temperatures above Tr (OsGM) and Tsr (RsGM). The 
results for both models are again indistinguishable for 
temperatures above Tr. Interestingly, mi is largely in- 
dependent not only of the system size, but also on the 
realization of the quenched disorder for the RsGM. 



C Other potential strengths 



We now turn to the question of the infiuence of the 
potential strength on the RsGM behavior. We have con- 
sidered two representative values, Vq = 0.2 and Vq — 5, 
i.e., five times smaller and larger, respectively, than the 
"canonical" value Vq = 1. The smallest value is close 
to that considered in Vq = 0.15, and we expect that 
our results will be comparable to theirs. As before, we 
begin by discussing the total roughness and the specific 
heat (see Fig. ^) . First of all, for all values of Vq there is 
a temperature above which the roughness value is inde- 
pendent of Vq and of the presence or absence of disorder. 
This means that our identification of this regime with the 
effective suppresion of any potential effect is indeed cor- 
rect: Different Vq leads only to different transition tem- 
peratures. Thus, for the OsGM we find r^«=o.2 ^ ^3 ± 1 
and T^'~^ = 19 ± 1, in agreement with the intuitive ex- 
pectation that larger potentials need higher temperatures 
to be suppressed. Aside from this, the general shape of 
the roughness curve is basically the same for the three 
values of Vq. The situation (for the OsGM) is the same as 
far as the specific heat is concerned: Larger (smaller) Vq 
leads to larger (smaller) anomalies, which are displaced to 
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higher (lower) T following the corresponding Tr. There- 
fore, we conclude that changing Vq does not introduce 
anything qualitatively new in the OsGM. 

The picture for the RsGM is substantially different: 
Modifying Vq does give rise to qualitatively new phe- 
nomena. Let us first look at the small Vq case. Figure 
^ shows that the roughness follows a straight line all the 
way down to T = (although we cannot exclude that 
there are nonlinear effects for T < 1 with our present 
resolution). This would suggest that there is no transi- 
tion in this case, very much like the results of Batrouni 
and Hwa Q . The upper panel of Fig. 0, where the height 
difference correlation function is depicted, confirms this 
interpretation, showing no dependence on temperature 
in the analyzed range; Fig. ||, for mi, agrees with this as 
well, in so far as the dependence of mi on temperature 
is well described by the high temperature expansion. In 
view of this, we can conclude that if there is a transition, 
it occurs at a temperature smaller than T ~ 1. 

Finally, let us consider the large Vg case, Vq = 5. The 
plot of the roughness in Fig. || exhibits a striking peak for 
T = 5, after which the roughness decreases until reaching 
the high temperature regime (marked by the correspond- 
ing OsGM with the same value for Vb) at T = 8 ± 1. To 
assess the relevance of this peak in the roughness, we 
repeated our simulations for L — QA and performed ad- 
ditional ones for L = 128. The results, collected in Fig. 

, show that the peak is a realization dependent feature. 
However, in this plot we also see that for temperatures 
below T — 5 the roughness is roughly independent of the 
system size, something which we did not observe when 
Vq — I (the lines for Vq = 1 are included in Fig. || again 
for comparison). Hence, even if the peak at T = 5 does 
not actually exist, that temperature does seem to sepa- 
rate two different regions. In addition, for Vq ~ 5 the 
specific heat has a (more smeared) maximum at about 
the same temperature as that of the roughness, although 
our data are much noisier and we cannot establish clearly 
the maximum temperature; the dependence, however, is 
manifestly non-monotonic. Figure ^supports our conclu- 
sion that T^1^~^ = 8 ± 1, whereas nothing special is seen 
as mi goes through T = 5, the roughness maximum. The 
most intriguing result is the one in Fig. ^ for the height 
difference correlation function: For T < 5, the scaled cor- 
relations decrease with temperature but, simultaneously, 
the correlation length increases, up to T = 5, when it 
increases beyond the system size. Above that tempera- 
ture, it follows the same evolution as the Vq = 1 case, 
finally reaching = 8 ± 1 . Whereas in this intermedi- 
ate temperature regime the height difference correlation 
functions are well described by squared logarithms. Fig. |^ 
immediately shows that the lowest temperature correla- 
tions can by no means be considered squared logarithms. 
This suggest the presence of a new phase transition at 
T* = 5 ± 1. We will discuss the possible nature of the 
low temperature phase and the existence of this transi- 
tion in the next section. 



D Other disorder strengths 

We now generalize the RsGM and let h^^\r) to be ran- 
domly chosen from a uniform distribution in [0, e], with 
< e < 2tt; e = 27r is the case studied in the preceding 
Section. As representative examples we have considered 
e = 0.8 and e = 0.2, closer to the RsGM and the OsGM, 
respectively. From the roughness dependence on tem- 
perature, shown in Fig. |l^, we see that for both values 
of e the dependence of the roughness on temperature is 
qualitatively similar to the OsGM, the results for the 
lower e value being practically the same as for the e = 
case. However, the case e = 0.8 is somewhat different: 
The low temperature region appears to consist of two 
straight lines, changing slope at a temperature around 
T' — 5, rather than a nonlinear dependence. By rea- 
soning as above, we identify Tg^-^ ~ Tr = 16 ± 1 and 
'^SR ^ = 12 ± 1, values which are confirmed by the en- 
ergy behavior (not shown), the height difference (Fig. [Tl| ) 
and slope (not shown) correlation functions, and by the 
dependence of mi on temperature (not shown). 

An interesting question arises from Fig. ^: There is no 
evidence about the squared logarithmic behavior found 
for the RsGM and, furthermore, the plots exhibit a fi- 
nite correlation length below Tsr, for both values of e. 
Another intriguing fact is the nonmonotonic dependence 
of the correlation function on temperature for e = 0.8: 
From the curve for T = 1, the scaled correlation func- 
tion decreases up to T = 5; upon further increasing the 
temperature, the evolution of the curves is very simi- 
lar to that of the OsGM. This might be connected with 
the change in slope in the roughness curve mentioned in 
the preceding paragraph (see Fig. |l0[ ), but we have not 
been able to draw a clearer connection. All this is clear 
evidence that the behavior of the RsGM is significantly 
dependent on the disorder strength. 



IV. DISCUSSION AND CONCLUSIONS 

Let us begin the discussion of the above results by an- 
alyzing our findings about the OsGM. Our simulations 
strongly support that Tr = 16 ± 1 for the OsGM on a 
lattice with Vq = 1. This is in contrast to the claims in 
||9|, |0) that Tr = 25 where a different way of definin- 
ing the transition, which assumes the validity of the RG 
approach, was used (see Sec. HI B and |4^, ^). Fur- 
ther, the result is also in contrast with the RG prediction 
itseff H, 1^, |l[ H, m, which in our units is Tr = 87r. 
However, we believe that the comparison with the linear 
theory for the EW high temperature phase has a much 
more physical character while keeping the basic RG ideas, 
and establishes beyond doubt that for the studied lattices 
the roughening temperature is Tfj = 16 ± 1 for Vq = 1. 
Another hint in favor of our claim is the finding of the 
(Schottky) anomaly in the specific heat, which should 
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appear below the transition temperature in view of what 
occurs for the XY and related models Finally, the 
fact that we obtain the same results for both the OsGM 
and the RsGM above Tr is clear evidence that the poten- 
tial is irrelevant (in the RG sense) in that regime and that 
we have indeed located the transition. Clearly, we can- 
not exclude the possibility that working on even larger 
lattices we would find the transition where the RG pre- 
dicts it, but the absence of any finite size effects even for 
L = 256 makes this possibility quite unlikely. Another 
possible reason for the discrepancy is the fact that our 
simulations are intrinsically discrete in space while RG 
theories for the OsGM are always applied to the contin- 
uum equation; again, very much larger lattices would re- 
move this objection and clarify the effects of discreteness. 
Aside from that, we have also found that increasing (de- 
creasing) Vq increases (decreases) the roughening temper- 
ature: In the cases we studied, we found T^°=°-^ = 13±1 
and T^'~^ = 19 ± 1, which is intuitively reasonable as 
larger potential barriers require larger temperatures for 
the surface to overcome them. On the other hand, RG 
calculations are perturbative in Vq, so one would expect 
better agreement to the RG prediction for Vq — 0.2, but 
in fact the agreement is worse in that case. 

Let us now turn to the RsGM. In the "canonical" 
case, Vb = 1, wc found a superroughening transition at 
TsR = 4 ± 1 = Tn/A, to be compared to RG predictions 
that it should occur at Tij/2. Below Tsr, we have ob- 
tained a In^ dependence of the height difference correla- 
tion function, in agreement with RG results. However, we 
have clearly shown that the superroughening transition 
temperature depends on Vb, confirming the earlier report 
by Batrouni and Hwa on the absence of the transition 
in Langevin dynamics simulations for Vq = 0.15, and 
later reports by Rieger [|j and Ruiz-Lorenzo who also 
observed this dependence in Monte Carlo simulations. In 
the opposite case, Vb = 5, wc find that T^^"^ = 8 ± 1, 
considerably higher than the Vb = 1 temperature. This 
disagrees with the RG predictions of a universal Tsr in- 
dependent of Vb . We believe that the agreement between 
our results and the previous ones |, |, indeed sup- 
ports a dependence of Tsr on Vb, whose explanation re- 
mains an open question as far as RG is concerned. 

A second, novel finding arises when considering our 
numerical results for Vb = 5, which strongly suggest 
the possibility of two different low temperature phases. 
In our comments in the preceding paragraph, we took 
Tg'j^^ — 8 ± 1 interpreting that the superroughening 
transition implies a change from a In^ behavior of the 
height difference correlation function to a In form (and 
the rest of EW features). However, the lack of size depen- 
dence of the roughness and the specific heat on tempera- 
ture below T* = 5 ± 1, with peaks absent for smaller val- 
ues of Vb, raises the possible existence of another phase 
transition. If one looks at the correlation functions in 
Fig. ^, it turns out that for temperatures below T* = 4 
the correlation length is finite, in agreement with the 



roughness independence on the system size. Whereas the 
range of correlations above T*, which we believe is infi- 
nite, could be a subject of debate as we only have studied 
sizes up to L = 256, our claim of finite correlation lengths 
below T* is difficult to dispute. Further evidence in this 
regard is shown in Fig. |l^, where curves for Vb = 5 at 
T = 1 are compared for two different system sizes. At 
this point, it is interesting to recall that in a previous 
paper it was found that at T = 0, the RsGM with 
Vb = 1 exhibits a finite correlation length of about 20 
lattice units (the reader may find Figs. 2 and 3 of Q 
illustrative). Having this in mind, it is not unreasonable 
to conjecture that there is a T* for the Vb = 1 case, which 
could be below T = 1 or close to 1 (the upper curve in 
Fig. m might already show a finite correlation length). We 
stress that this phase has not been pre viously reported in 
works at T = 0: Thus, Zeng et al. studied a discrete 
model but, being different from the Gaussian, their low 
temperature results can not be compared to ours, and the 
results of Rieger et al. [|l5|, |l^ do not allow to conclude 
anything in this respect. Intuitively, one can expect a 
finite correlation length phase at low temperatures and 
large Vb; in the limit Vb — > cxd, the surface follows the 
disorder (i.e., h{r) = h{r)^°'> + 2n(r)7r everywhere), but 
the gradient term smoothens out the lack of correlations 
of h^^^'lr), the competition of these two effects yielding a 
finite correlation length. In a loose sense, this could be 
interpreted as Anderson localization taking over the cou- 
pling between neighboring sites with increasing Vb . This 
picture is confirmed by simulations for Vb = 25 (Fig. 
For such a large value of Vq the correlation length is only 
1 lattice unit, i.e., correlations reach only nearest neigh- 
bors. Clearly, the data presented here is not conclusive, 
but the conjecture that there are two transitions whose 
critical temperatures depend on Vb is not unreasonable 
and deserves further consideration. 

We have also found that the transition temperature 
and the correlation functions depend on the disorder 
strength. This is not unexpected, in so far as the change 
in the disorder distribution interacts with the periodicity 
of the sine potential, and therefore it is clear that when 
e — 2tt, i.e., in the standard RsGM, we are in a special 
case. In this respect, our results appear to indicate that 
the RsGM (with e = 27r) is a very specific model, and that 
its behavior at low temperatures might not be represen- 
tative of what one would find in an actual experiment, 
where the disorder can not be so precisely controlled. 
Another conclusion we may draw from our work is that 
there might be two classes of behavior at low tempera- 
tures for small and large e: Small e models would behave 
very similarly to the OsGM, whereas large values of e 
would give rise to a more complex phenomenology with, 
e.g., nonmonotonic behavior of the correlation functions. 

As a final conclusion, we remark that the most rele- 
vant result of the present work is the determination of the 
transition temperature from the high temperature phase 
to the low temperature phase (or phases) of both the 
OsGM and the RsGM. This poses a number of questions 
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to be addressed either with greater numerical capabiUties 
or with new analytical tools. We believe that the complex 
phase diagram of the RsGM is being partially unveiled, 
as our research supports previous findings such as the 
non universality of the transition temperature. Signifi- 
cantly, once we know where to look for the low tempera- 
ture phase of the RsGM, we can investigate the nature of 
that phase (or phases). It is often stated that the RsGM 
is "glassy" in this regime; however, this assertion has 
never been really proven nor fully detailed and, further- 
more, if the RG picture is qualitatively correct, it has to 
be recalled that it is a replica-symmetric theory, which 
implies that the superrough phase would not be glassy in 
the replica sense . We have obtained preliminary ev- 
idence that there are long-lived metastable states in the 
low temperature phase of the RsGM but in view 
of our present results and their non-universality we will 
examine this question more carefully in future work. In- 
vestigation of the dynamics of the RsGM will also be 
important; we recall that Batrouni and Hwa Q did find 
evidence for a superroughening transition in the dynam- 
ics of the model, and hence it would be worth checking 
their results for larger values of Vq. We hope that those 
analyses, along with measurements of nonlinear suscepti- 
bilities and of relaxation dynamics, will shed light on this 
difficult problem. Work along these lines is in progress. 
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FIG. 1: Mean energy for the OsGM and the RsGM vs temperature. Symbols as indicated in the plot. The straight 

lines correspond to the low temperature prediction (lower line) and the high temperature prediction (upper line). 
Inset: Specific heat vs temperature; symbols as in the main plot. Error bars are smaller than the symbol size. 




FIG. 2: Total roughness for the OsGM and the RsGM vs temperature. Symbols as indicated in the plot. The straight 
lines correspond to the high temperature prediction for the two sizes considered. Inset: Roughness susceptibility vs 
temperature; symbols as in the main plot. Examples for error bars in different regions are shown. 
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FIG. 3: Height difference correlation functions (scaled by T) vs Inr, r being distance, for the OsGM (upper panel) 
and disordered (lower panel) sine-Gordon model. Temperatures are indicated at the right side of the plots. The curve 
marked with temperature 16 in the plot for the RsGM is the correlation function for the OsGM at that temperature, 
showing clearly that both overlap. 
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FIG. 4: Slope correlation functions (scaled by T) vs Inr, r being the distance, for the OsGM (upper panel) and 
disordered (lower panel) sine-Gordon model. Temperatures are indicated at the right side of the plots. The curve 
marked with temperature 16 in the plot for the RsGM is the correlation function for the OsGM at that temperature, 
showing clearly that both overlap. 
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FIG. 5: Comparison of the averages of the cosine term, mi, of the OsGM and the RsGM at different temperatures. 
Symbols and hnes as indicated in the plot; lines correspond to theoretical approximations up to order j3 and 
whereas symbols are numerical results. Results are independent of the system size and of the realization of the 
disorder; error bars for thermal averages are smaller than the symbol size. The inset is an enlargement of the low mi 
region of the plot. 
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FIG. 6: Comparison of the roughness (upper panel) and specific heat (lower panel) for the OsGM and the RsGM for 
different values of Vq. The inset in the upper panel shows a larger range of temperatures for the roughness dependence. 
Symbols as indicated in the plot. All the results have been obtained for L = 64. Note especially the peak in the 
specific heat for the RsGM with Vq = 5. 
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FIG. 7: Comparison of the height difference correlation functions for the RsGM with Vo, Vo = 0.2 (upper panel) and 
Vq = 5 (lower panel). Temperatures as indicated in the plots. 




FIG. 8: Comparison of the averages of the cosine term, mi, of the OsGM and the RsGM for different values of Vq. 
Symbols as indicated in the plot, and lines correspond to the theoretical approximations for high temperature. 
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FIG. 9: Comparison of the roughness OsGM and the RsGM for Vq = 1> 5 and sizes L = 64, 128. Symbols as indicated 
in the plot. 
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FIG. 10: Comparison of the roughness for the OsGM, the RsGM, and the two intermediate versions of the RsGM 
with disorder strength e = 0.2 and e = 0.8 (see text for the corresponding definition). The size is L = 64 in all four 
cases, and symbols are as indicated in the plot. 
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FIG. 11: Height difference correlation functions (scaled by T) vs Inr, r being distance, for the e = 0.2 (upper panel) 
and f = 0.8 (lower panel) RsGM. Temperatures arc indicated at the right side of the plots. The curve marked with 
temperature 16 in the plot for the RsGM is the correlation function for the OsGM at that temperature, showing 
clearly that both overlap. 
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FIG. 12: Height difference correlation functions (scaled by T) vs Inr, r being distance, for different values of Vq and 
L, as indicated in the plot. In all cases, T = 1. 



